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We investigate tiie phase portrait of the (H-l)-dimensional massless two-flavored NJL2 model 
containing a quark number chemical potential /x and an isospin chemical potential in the limit 
of a large number of colors — >■ 00. Particular attention is paid to the question to what extent 
the inclusion of an isospin asymmetry affects chiral condensates to have a spatial inhomogeneity 
, in the form of the so-called chiral density waves (CDW) (chiral spirals). It is shown that at zero 

• temperature and comparatively small values of ^, i.e. at < /Xc ~ O.68M0 (Mo is the dynamical 

' quark mass in the vacuum) only the homogeneous charged pion condensation phase is realized for 

^sj , arbitrary nonzero values of /jj. Contrary to this, for large values of /i > /^c, two CDW phases appear 

I ■ in the {fii , ^)-phase diagram of the model. In the first phase, CDWs are clockwise twisted chiral 

^ ' spirals and in the second phase they are counterclockwise. The influence of nonzero temperature on 

^~~> , the formation of the CDW phases is also investigated. 

\o ■ 

PACS numbers: 12.39.Ki, 12.38.Mh, 21.65.Qr 
Ph' I. INTRODUCTION 

1^ ' During the last decade, much attention has been attracted to the investigation of the QCD phase diagram in terms 
of quark number as well as isospin chemical potentials. First of all, this is motivated by heavy-ion collision experiments 
where dense baryonic matter has an evident isospin asymmetry, i.e. different neutron and proton contents of initial 
ions. Moreover, the dense hadronic/quark matter inside compact stars is also expected to be isotopically asymmetric. 
Generally speaking, it is understood that one of the important challenges for QCD is to describe the dense and hot 
baryonic matter in different physical situations. However, in the above mentioned realistic situations the quark density 
is rather small, and weak coupling QCD analysis is not applicable. So, different nonperturbative methods or effective 
theories such as chiral effective Lagrangians and especially Nambu - Jona-Lasinio (NJL) type models (l| are usually 
employed for the consideration of the properties of dense and hot baryonic matter under the conditions of heavy-ion 
experiments or in the compact stars interior, i.e. in the presence of external factors such as temperature, chemical 
potentials, magnetic field, finite size effects, etc. (seeje.g., and references therein). In particular, phenomena of 
■ dense quark matter like color superconductivity [2, Q-Q , as well as charged pion condensation P-flTj were investigated 
• • ' in the framework of these QCD-like effective models. 
. 5^ . It should be noted that an effective description of QCD in terms of NJL models, i.e. through an employment 
I of four-fermionic theories in (3-|-l)-dimensional spacetime, is usually valid only at comparatively low energies and 
;h • densities. At the same time, (l-l-l)-dimensional Gross-Neveu (GN) type models p^. [T3| are valid also at high energy 
and density, and due to their properties of renormalizability, asymptotic freedom and spontaneous chiral symmetry 
breaking, can also be used for a reasonable qualitative modeling of QCD even at finite temperature and hadron density 
p^ - [l8| . Due to the relative simplicity of GN models in the leading order of the large TVc-expansion {Nc is the number 
of colored quarks), their use is convenient for the application of nonperturbative methods in quantum field theory [1^ . 
Moreover, it is worth noting that it is in the leading order of the large iVc-expansion that the well known no-go theorem 
of Mermin-Wagner-Coleman plj , apparently forbidding the spontaneous breaking of continuous symmetries in the 
(l-l-l)-dimensional models, becomes invalid |15l - [l8| . (It means that in the large Nc limit quantum fluctuations, which 
would otherwise destroy a long-range order corresponding to a spontaneous symmetry breaking, are suppressed by 
l/Nc factors.) Note also that GN type models are quite suitable for the description of physics in quasi one-dimensional 
condensed matter systems such as polyacetylene |20| . 

Thus, such phenomena of dense QCD as color superconductivity, where the color group is broken spontaneously, 
and charged pion condensation, where spontaneous breaking of the continuous isospin s ym metry takes place, might be 
modeled in terms of renormalizable (l-l-l)-dimensional GN type models (see, e.g., [l7j and [23 - |23 |. respectively). 

In our previous papers [2342^ the phase diagram of a (l-l-l)-dimensional S'C/l(2) x SUr{2) symmetric NJL model 
with two massless or massive quark flavors was investigated in the leading order of the l/7Vc-expansion and in the 
presence of the quark number- as well as isospin chemical potentials. There we have considered the case of order 
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In this paper we shall use the notation NJL model for theories with four-fermionic interactions also for (l+l)-dimensional models with 
a continuous chiral symmetry group instead of "chiral 2D GN model" due to the fact that the chiral structure of the Lagrangian is 
indeed closely related to the (3-|-l)-dimensional NJL model. 
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parameters (condensates) that are homogeneous, i.e. independent of the space coordinate. The situation corresponds 
to the conserved Lorentz and spatial translational invariance and is adequate to physical systems in vacuum, i.e. 
at zero chemical potentials. In dense baryonic matter, i.e. at nonzero quark number chemical potential, there 
might appear new phases with a spatially inhomogeneous chiral condensate which destroys both chiral and spatial 
translational invariance of the syste m ( see the relevant discussions made in the framework of both (l+l)-dimensional 
[3 [2^ HB] and (3+l)-dimensional j27l - l3ll | models). Thus, in this paper and in contrast to [23 - l24 |. we consider the 
phase portrait of the above mentioned massless SUl (2) x SUr{2) symmetric NJL model with two chemical potentials in 
the leading order of the 1 /TVc-expansion taking into account the possibility that the chiral condensate might become 
inhomogeneous and take the form of a (dual) chiral density wave (CDW). In this case the scalar quark-antiquark 
condensate, {qq), and the pseudoscalar condensate of neutral 7r° mesons, {q'y^T^q), form a chiral spiral, i.e. 

{qq) ^ cos26x, {qi^T^q) ^ sin26x, 

where T3 is the isospin Pauli matrix, x is the space coordinate, and 6 is a wave vector which has to be determined 
dynamically through the thermodynamic potential. It is necessary to point out that the inhomogeneous CDW- 
condensate is relevant to dense quark matter (32j and the chiral magnetic effect [s^. Both phenomena probably 
might be observed in heavy-ion collision experiments, where isotopic asymmetry is an inevitable property. So, we 
believe that the investigation of CDW condensates in the framework of the two-dimensional NJL model with isospin 
chemical potential could shed some light on physics of heavy-ion collisions. It should be noted, however, that our 
investigation of a condensate inhomogeneity in the form of the chiral density wave is only a first step. There may exist, 
at least at zero isospin chemical potential, other more preferabl e sp atially non-uniform ground state configurations of 
the chiral condensate like, for instance, chiral crystals [2^, [2^ [2^) (in the last case only the scalar {qq) condensate 
is an inhomogeneous quantity, and other condensates are homogeneous ones), but, in general, they are much harder 
to deal with. For technical reasons, in studying CDW configurations we do not take into account a nonzero bare 
(current) quark mass, although recently some efforts to get rid of this assumption have been made (30| . 

The paper is organized as follows. In Section II, we derive in the leading order of the large iVc-expansion the 
expression for the thermodynamic potential of the SUl{2) x SUii{2) symmetric massless NJL2 model with quark 
number chemical potential fi and isospin chemical potential at zero temperature. Here we also consider the 
possibility of a spatial inhomogeneity for the chiral condensates in the form of the so-called chiral density waves. 
First, the phase portrait of the model is discussed in the simple case of spatially homogeneous condensates in Section 
III, and then, in Section IV, the phase structure of the model in terms of /i and fij and at zero temperature is 
investigated for CDW inhomogeneous phases. The influence of nonzero temperature on the formation of CDW phases 
is considered in Section V. Finally, Section VI presents some concluding remarks. 



II. THE MODEL AND ITS EFFECTIVE ACTION 



We consider a (l+l)-dimensional NJL-type model with two massless quark flavors {u and d quarks) to mimic 
properties of real dense quark matter. Its Lagrangian has the form 

{qqf + iqil'rqf], (1) 

where our choice for the gamma matrices in (l+l)-dimensions is as follows: 7° = 0-2,7^ = io'i,7'^ = 7°7^ = cs, and 
the quark field q{x) = qia{x) is a flavor doublet (i = 1,2 or i — u, d) with corresponding Pauli matrices (fc = 1, 2, 3) 
and color iVc-plet (a — 1, ...,Nc), as well as a two-component Dirac spinor (the summation in ([1]) over flavor, color, 
and spinor indices is implied). The quark number chemical potential in ([T]) is responsible for the nonzero baryonic 
density of quark matter, whereas the isospin chemical potential /x/ is taken into account in order to study asymmetric 
quark matter at nonzero isospin densities (in this case the densities of u and d quarks are different). Evidently, the 
model (dJ is a generalization of the original (l-l-l)-dimensional Gross-Neveu model [l^ with a single quark to the 
case of two quark flavors and additional chemical potentials. As a result, we have for our model a more complicated 
chiral symmetry group. Indeed, at /i/ = apart from the global color SU(A'c) symmetry, the Lagrangian ([T]) is 
invariant under transformations of the chiral SUl{2) x SUr{2) group. However, at fij this symmetry is reduced 
to Ui^l{1) X Ui^r{1), where I3 — T3/2 is the third component of the isospin operator (as usual the subscripts L,R 
mean that the corresponding group acts only on the left, right handed spinors, respectively). Evidently, this symmetry 
can also be presented as Ui^{l) x UAisil), where [//^(l), [7^/3(1) denote the isospin and the axial isospin subgroups 
respectively. Quarks are transformed under these subgroups as g — >■ exp(iar3)q and q — t- exp(ia7''^r3)g, respectively q 
Notice that Lagrangian ([1} is parity invariant. 



£ = q 



7Z 



^ Recall that exp(iaT3) = cos q + its sin a, exp(io7^T3) = cos o + i7^T3 sin o. 
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The linearized version of the Lagrangian ([T]), which contains composite bosonic fields a{x) and TTa{x) (a = 1,2,3), 
has the following form: 



C 



Yidp + ^7" + J^T37° - CT - ij^TTaTa 



El 

4G 



(7 +TT„ 



(2) 



where v = /^//2. Evidently, the Lagrangian ([2]) is equivalent to ([T]), which simply follows from the use of the following 
constraint equations for the bosonic fields 

G 



a{x) 



'2— (gg); ^aix) 



Q 



(3) 



Furthermore, it is clear from ([3]) and footnote [5] that the bosonic fields transform under the isospin Ui^{l) and axial 
isospin C/4/3(l) subgroups in the following manner: 



[7/3(1) : a ^ a\ tts ^ tts; tti — > cos(2a)7ri + sin(2Q!)7r2; 1^2 — ^ cos(2a)7r2 — sin(2a)7ri, 
[7^/3(1): TTi — > vTi; 7r2 — > 7r2; — > cos(2a)CT + sin(2Q)7r3; tts cos(2a)7r3 — sin(2a)(T. 



(4) 



Starting from Lagrangian ([2]), one obtains in the leading order of the large A'^c-expansion (i.e. in the one-fermion loop 
approximation) the following path integral expression for the effective action iSeff (u, tTq) of the bosonic a(x) and 7ra(a;) 
fields: 



exp(i5efi(CT,7ra)) = iV' / [(i(j][dg] exp(i / td^x], 



where 



5cff(o-, vTa) = -TVc / da; 



1 d^a; 
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+ 5, 



cff , 



(5) 



and N' is a normalization constant. The quark contribution to the effective action, i.e. the term 5cff in ([5]), is given 
by: 



exp(. 



(i5eff) -iV' /"[rfg][dg]cxp(i / {g[7''iap + Ai7° + i/r37"-fT- 17571, Tajgjd'a;). 



(6) 



The ground state expectation values (cr(a;)) and (7ra(x)), of the composite bosonic fields are determined by the saddle 
point equations. 



Sa{x) 



0, 



S-Kaix) 



0, 



(7) 



where a — 1,2,3. In vacuum, i.e. in the state corresponding to an empty space with zero particle density and zero 
values of the chemical potentials /j, and /i/, the quantities (c^x)) and (iTaix)) do not depend on space coordinates. 
However, in a dense medium, when ^ ^ 0, /i/ ^ 0, the ground state expectation values of bosonic fields might have a 
nontrivial dependence on x. In particular, in this paper we will use the following ansatz: 



{a{x)) = M cos{2bx), (713(2;)) = M sin(26a;), (7ri(a;)) 



(Mx)) = 0, 



(8) 



where M, b, and A are constant quantities. In fact, they are coordinates of the global minimum point of the thermo- 
dynamic potential (TDP) n{M,b,A). In the leading order of the large A^c-expansion it is defined by the following 
expression: 



d'^xn{M, b,A) 



1 



ScsWix), TTaix)} 



a{x)—{a{x)) ,7Ta{x)—{7ra(x)) ' 



which gives 



where 



i / d^xn{M,b,A) ^ i d^x 



A,f 2 I A2 1 

- In ( / idq] [dq] exp ( i / d'xqVq 



AG 



Nr. 



(9) 



(10) 



V = Y'^dp + /i7° + z^r37" - M exp(2i7'^r36x) - i^i^TiA. 



(11) 



Here and in the following we will use a rather conventional notation "global" minimum in the sense that among all our numerically 
found local minima the thermodynamical potential takes in their case the lowest value. This does not exclude the possibility that there 
exist other inhomogeneous condensates, different from JSj, which lead to ground states with even lower values of the TDP. 
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To proceed, let us introduce the new quark fields, q„, = exp ( 17^x3 6a; )(7 and = 9 exp(i7^T36a;), such that 

qVq = q^ [^idp + fl-f"^ + (b + iy)T3j'^ -M- i-f^nA] q^ = q^Dq^, (12) 

where instead of the x— dependent Dirac operator (jlll) a new x— independent operator appears 

D = Y'^dp + + + i^)t37" -M - i7'ViA. (13) 

Since this transformation of quark fields does not change the path integral measure in ([T0|) 0, expression (fTO|) for the 
thermodynamic potential is easily transformed into the following one: 

, AX Af2 + A2 .Tr^f^lnZ) 

= + iTr./ j j^\n[^ + m' + {b + v)T,f -M- 17^ An) , (14) 

where the Tr-operation Tisfx stands for the trace in spinor- (s), flavor- (/) as well as two-dimensional coordinate- 
{x) spaces, respectively, and Tr^/ is the respective trace without x— space. Since the thermodynamic potential (|14p is 
formally equal to the TDP (9) of paper [23] when one performs the replacement v ^ h + v, one can further use the 
corresponding techniques and obtains 

n{M, 6, A) = + V (^ { [(^0 + (^a)'] [bo + m)' - {E-^f\ }, (15) 

where 



- + =i;±(6-hz^), E = ^pI + M^. (16) 

The argument of the ln(a;)-function in ()15p is proportional to the inverse quark propagator in the energy-momentum 
space representation. Hence, its zeros are the poles of the quark propagator. So, using ([15]) one can find the dispersion 
laws for quasiparticles, i.e. the momentum dependence of the quark (pou, Pod) and antiquark {poa, Pnd) energies, in a 
medium (the full expression of the quark propagator matrix is presented in Appendix B of paper ^23] ) : 

P()u = E^-fi, pM = E^-fi, pou = -{EX+ fJ.), Paa^ -{Ea+ f^)- (17) 

It is clear that expression (jisp is symmetric with respect to the transformations fi — )• — /i and {b + v) — > —(b + v) 
respectively. Thus, without loss of generality, it is sufficient to adopt the restrictions /i > and {b + i') > 0. Under 
these conditions, upon integrating in ([TSl) over pq (see in [l^l for similar integrals), one obtains for the TDP of the 
system at zero temperature the following expression: 

+ EX)e{ii - E+) + E^)e{fi -e^)}, (i8) 

where 9{x) is the Heaviside theta- function. 

III. HOMOGENEOUS CHIRAL CONDENSATE, fe = 
A. The case with flavor symmetry, /i/ = 

First of all let us consider the vacuum case, i.e. when /j, = 0, /u/ = and temperature is zero. Since in QCD parity 
is not broken in the vacuum, it is necessary in all QCD-motivated theories to adopt the same requirement. In the 
framework of our model this means that one should take A = if /x = 0, /x/ = 0. Assuming homogeneity of the chiral 
condensate (6 = 0), we then obtain from (|18l) the following expression for the effective potential of the initial NJL2 
model in vacuum (A = 0, /x = 0, /u/ = 0) 0: 

Vo{M) = — --j^ dpr^pj + M2. (19) 



* This nontrivial fact follows from the investigations by Fujikawa i34il . who established that a chiral transformation of spinor fields changes 

the path integral measure only in the case when there is an interaction between spinor and gauge fields. 
^ In vacuum, the thermodynamic potential is usually called effective potential. 
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Formally, the effective potential (IT9|) is an ultraviolet (UV) divergent quantity. To renormalize Vo{M), i.e. to obtain a 
finite expression for it, we first need to regularize the integral in the right hand side of (jl9l) by cutting off its integration 
region, pi < A. Second, we suppose that the bare coupling constant G in ([T9|) depends on the cutoff parameter A 
(G = G(A)) in such a way that in the limit A — > oo one obtains a finite expression. To construct the function G(A), 
let us suppose that the stationarity equation dVo{M)/dM = has a nontrivial solution Mq. Then it is easy to obtain 
from this equation the following expression for the bare coupling constant G(A): 

12'-^ 




dpi , . --In . (20) 



2G(A) nJo "'VW 

Now, using (|20p in the regularized expression (flQl) and adding an unessential constant A^/tt, one can find at A — ^ cx) 

f M2 



M 2 f I 

|4G(A)-^X ^^^a/^ 



A2' 


M2 












m 




- 1 


.J 


27r 





(21) 



Since Mq might be considered as a free model parameter, it follows from ([20)1 and (|2ip that the renormalization 
procedure of the NJL2 model is accompanied by the dimensional transmutation phenomenon. Indeed, in the initial 
unrenormalized expression (jl9p for Vo(M) the dimensionless coupling constant G is present, whereas after renormal- 
ization the effective potential (|2T|) is characterized by a dimensional free model parameter Mq. Moreover, as it is clear 
from (|2ip . the global minimum point of the effective potential Vq{M) lies just at the point M — Mq, so in vacuum 
the chiral SUl{2) x SUr{2) symmetry of the NJL2 model (1) is always spontaneously broken and the quantity Mq 
might be treated as dynamical quark mass (in vacuum) . 

A detailed information about the phase structure of the NJL2 model (1) at ^ 7^ 0,^/ = can be found, e.g., in (141] . 
So, at /i > Mq/\/2 and fj,i — there is a massless chirally symmetric phase with nonzero baryon density. However, at 
fi < Mq/^/2 and ni — chiral symmetry is spontaneously broken down and quarks acquire a mass Mq. In this phase 
baryon density is equal to zero. 



B. Phase structure in the general case: /i 7^ 0, /i/ 7^ 



To find the phase portrait of the NJL2 model (1) in the case of a homogeneous chiral condensate but for arbitrary 
values of chemical potentials and at zero temperature, one should start from the expression ^TE\i with 6 = 0. (Note, 
that at /i/ 7^ the condensation of charged pions might occur, so we need to take into account a nonzero value of 
A.) Obviously, this expression is again UV-divergent, so first of all it is necessary to regularize it. Using, as the most 
simple regularization, a A-cutoff in the one-dimensional momentum space, we have: 



f}rog(Af,6 = 0,A) 



M' + 
4G 



dpi f 



F+ 



-■} 



^{(/i - £+)9{f^ - £+) + (a* - £^)0(/i -£Z)}, 



(22) 



where £^ denotes the quantity ([T6l) at & = 0. Due to the presence of (^-functions, the second integral in ([22l) has 
a finite integration region, i.e. it is a proper integral not needed to be regularized. To obtain a finite (renormalized) 
expression ri(M, A) for the thermodynamic potential, one should again perform in the replacement G — G(A), 
the last quantity being given in ((20)) . and then let A tend to infinity (compare with ([2T|)). i.e. 



n{M, A) = lim <^ n.^JM, 6 = 0, A) 

A— )-oo 



G^G(A) 



A2 

TT 



(23) 



Using the definition of the effective potential in vacuum (see (pij) '). it is easy to obtain the following renormalization 
invariant expression of the TDP (|23)) 



dpi (r.+ 



£; 



pI 



n{M, A)= Vq{^M^ + A2) - 

- ^{(m - £mt^ - £l) + (m - £^)e{^^ £^)}, 



M2 + A2 



(24) 



where the function Vq{x) is defined in ([21]) . Moreover, the second integral in (|24p is proper (see also the corresponding 
remark just after ([H])), whereas the first integral is convergent and defined as 



dpi 




£\ + f A 



M2 



(25) 
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FIG. 1. The (/i, v) phase portrait of the model considered 
at T = and v > Q in the case of spatially homogeneous 
condensates. Here v = and Mo is the quark mass 
in the vacuum. Number 1 denotes the symmetric phase 
with massless quarks, number 2 the normal quark matter 
phase with massive quarks, and PC denotes the charged 
pion condensed phase. The point a is the lowest point of 
the phase 2: ~ O.68M0, « O.6M0. 



FIG. 2. The plot of QP^^'' O vs M, h at /i/ 
Mo, A = 0. 



0,M 



Thus, in the case of a homogeneous chiral condensate, the TDP is given by ([23|) - p4|) and the corresponding phase 
structure, following from it, is depicted in Fig. 1 (for a more detailed investigation of this TDP see paper p^l. 
There, the phases denoted by 1 and 2 correspond to the global minimum point (GMP) of the form (M = 0, A = 0) 
and {M 7^ 0, A = 0), correspondingly. In the pion condensed phase (PC) the GMP of the TDP has the form 
{M = 0, A = Mq), i.e. in this phase the isospin symmetry Ui,^{l) is broken spontaneously 0. It is easy to see that 
throughout the PC phase the quark number density is equal to zero, whereas the isospin density n/ — —dfl/diij is 
equal to v/tt. 



IV. INHOMOGENEOUS CHIRAL CONDENSATE, 6/0 



To obtain the phase portrait of the initial NJL2 model in this case (temperature T is zero), let us start from the most 
general expression for the TDP (|18l) . As previously, let us first use the most simple momentum cutoff regularization 
of this quantity. 



acg(M,6,A) 



4G 
dpi 



Jo I 



I E+)e{^, e+) + {^,- E-^)e{t, -ed], 



(26) 



where the expressions for E^ are presented in (|16p. The corresponding renormalized expression for the TDP is again 
defined by (compare with ([H])) 

n{M,b,A)= lim \ n,cg{M,b,A) 



G-^G(A) TT 



(27) 



where G{A) is given in d^H]) . and reads 

n{M, b, A)= Vb( + A2) - lim I / — 

Jo 



ei + e^-2^pI + ap + a? 
^{(m - EX)6{pi - E+) + (m - E-^)e{^i -E-^)]. 



(28) 



Note that our numerical investigations show that the TDP 1124 l l has no local minima of the form (M / 0, A ^ 0), i.e. simultaneous 
dynamical quark mass generation and charged pion condensation are incompatible in the framework of the NJL2 model (1) at b = 0. 
The same is valid for the simple two-flavored NJL4 model in the mean- field approximation [Tol |. However, it is not excluded that there 
might be realized a mixed phase with both nonzero gaps, M 7^ and A 7^ 0, in models with a more complicated four-fermion structure. 
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(Evidently, at & = this expression coincides with the TDP r2(M, A) (|24|) .1 The global minimum point of the function 
r2(M, b, A) vs variables M, b and A should render the phase structure of the model. However, two circumstances 
prevent us from considering this quantity as a genuine physical thermodynamic potential of the system. The first is 
that the function is not bounded from below with respect to the variable b. Secondly, it is intuitively clear that 
at M = the genuine thermodynamic potential should not depend on the variable b, because no observable quantity 
may depend on a wave vector if the amplitude of the corresponding oscillations (wave) is zero. However, the TDP 
defined by ^ at M = (see also in ^) 



n{M = 0,6, A) = Vb(A) 



(6- 



9{p-A) 



A^ In 



(29) 



retains an unphysical dependence on b. Clearly, the two above mentioned unphysical properties of the TDP 



due to the term 



in (P^ . Hence, the subtraction of this term from the TDP (|28p brings us to the quantity, 



which might serve as a physically acceptable thermodynamic potential of the system. 



nP^y''{M,b,A) ^n{M,b,A)■ 



ib+ly)' 



(30) 



(We also add in the expression pop a 5- independent term, -^^/tt, in order to reproduce at 6 = the TDP 
corresponding to a spatially homogeneous chiral condensate.) The reason why the expression for the TDP (f30|) 
does not follow straightforwardly from the unrcnormalized TDP expression (|18l) lies in the usage of the symmetric 
momentum cutoff regularized TDP (j26p . This means that for each energy the integration in the first (regu- 
larized) integral of (j26p is performed over the same momentum interval < pi < A. Correspondingly, in this 



case there is an asymmetry in values of energies E^, which contribute to flrcg{M,b, A). Indeed, if pi < A, then 



E 



< 



lvP±{b + v)') + A^, i.e. for different quasiparticles there are allowed different regions of their 
energy values. However, as discussed in the recent papers [28l . [2^ a more adequate regularization scheme in the 
case of spatially inhomogeneous phases is that one, where there is an energy constraint which is the same for all 
quasiparticles. So, dealing with spatial inhomogeneity, one can use, e.g ., the Schwinger proper-time regularization, 
dimensional regularization etc. In particular, in the recent paper [Slj the symmetric energy cutoff regularization 
scheme was proposed in considering the behavior of chiral density waves in the presence of an external magnetic field 
in the framework of a four-dimensional Nambu-Jona-Lasinio model. There, for each quasiparticle the same (finite) 
interval of their energy values was allowed to contribute to the regularized thermodynamic potential. As a result, a 
physically relevant renormalized TDP without the above-mentioned shortcomings was obtained. 

In this paper the slightly modified energy cutoff regularization scheme of is adopted. Namely, we require that 
only energies with momenta pi, constrained by the relations E^{M = , A — 0) — pi ± {b + v) < A, contribute to the 
regularized thermodynamic potential. This means that the term with energy E'^ {E~^) should be integrated in the 
regularized expression for TDP over the interval < pi < A — (6-|-J^) (0 < pi < A-|- {b + v)). (A similar regularization 
was used in studying the CDW phase in a two-dimensional NJL model without isospin chemical potential [13] ■) 
Consequently, we have the following regularized expression for the TDP 



a-cg(A/, 6, A) 



M2 



4G 



dpi Ej 



A+i> 



dpi E. 



Jo L 



E+)^{^l 



E+) + {y. - 



E^)e{^x^E-^)], 



(31) 



where v — {b + v). Replacing in this formula G by G(A) from (|20)) and adding an unessential constant (A^ — v'^)/-k, 
we obtain a physically "improved" renormalized expression Vl{M, b, A) for the TDP (|18l) when A — > oo, which differs 
from the expression n{M, b, A) in ([28| . Comparing (|26|) and pip one can easily find that 



n{M, b, A) 



rj(M, 6, A) = lim 

A— >cxo 



1 



A-P 



El 



1 



dpi E. 



{b- 



(32) 



(To obtain the last expression in p2p one should take into account that at A — > cx) the pi-values in both integrals are 
much greater than M, A, 6, y, pij. In this case it is possible to expand the quantities E^ into power series of pi and 
then to integrate each term.) Comparing (15(71) and ([5^ . we see that r2(M, &, A) = ^^^'^^{M ,b, A), i.e. there exists 
a regularization scheme Q, which in the case of the inhomogeneous chiral condensate ([5]) brings us to a physically 



Moreover, we expect that any regularization scheme, in which there is a constraint on the energy values common for all quasiparticles, 
should provide us with TDP f7P'^y''(M, 6, A) 13011 . Among these regularizations are dimensional and analytical ones, Pauli-Villars and 
Schwinger prope-time regularizations, as well as the above mentioned symmetric energy cutoff regularizati on l3ll. In particular, the 
proper-time regularization is often used in studying inhomogeneous phases in the framework of NJL models [281 . 12911 and does not lead 
to any unphysical effects, etc. 
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FIG. 3. The CDW amplitude M{fi) and its wave vector 
boifJ.) as functions of /x in the case of zero isospin chemical 
potential. Here Hc = Ha ~ O.68M0. 




FIG. 5. Section of the plot of QF^^"" (33) vs 6 at = Mo 
along the axis b passing through the point of minimum. 




FIG. 4. The plot of QP'^^'' (33) vs M, 6 at = Mc 
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FIG. 6. Section of the plot of QF^^^ (33) vs M at = Mo 
along the axis M passing through the point of minimum. 



acceptable TDP f]Pi>y''(M, &, A) Notice also that if 6 = then f7(M,6, A) is equal to 17P'^y"(M, 6, A). Hence, 

in the case of homogeneous chiral condensates the two above considered regularization schemes are equivalent. In 
contrast, in the inhomogeneous case the thermodynamic potentials Vl{M, b, A) and ^^^^^{M, b, A) differ by terms, 
containing the dynamical quantity b. As a result, the regularizations are not equivalent. However, since the symmetric 
momentum cutoff regularization (SMCR) is easier to handle, it is possible to perform all calculations in the framework 

of the SMCR scheme and then simply correct the obtained TDP 17(M, &, A) by the terms ^^^^ - ^ l^^^ (1501) ]. 
instead of using from the beginning one of the physically acceptable regularizations bringing us directly to the TDP 
np^y<^{M, b,A). 

To illustrate the fact that the TDP fiPi^y" (M, 6, A) is bounded from below as a function of the variable b, we plot 
in Fig. 2 this thermodynamic potential vs M, b in the particular case /z/ = 0, A = 0, = Mo. 



A. Particular case: = 0, /i 7^ 

Recall that the CDW inhomogeneous phase was established earlier in the NJL2 model with ?7l(1) x [/^(l) chiral 
symmetry for all /i > at rather low temperatures [itI. [18|. In contrast, in this paper we are going to study chiral 
density waves in the NJL2 model with a continuous chiral SUl{2) x SUr{2) symmetry. In the present section we 
consider the case of T = 0. It is well known that at /i/ = the charged pion condensation phenomenon is forbidden 
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(see, e.g., in [IJl), so without loss of gen erality one may suppose that A = in (|28l) . Then the TDP r2(M, &, A = 0) can 
be easily evaluated analytically (see [24|) and the physical thermodynamic potential fl^^^'^^lM, b) = fl^^^^lM, b,A = 0) 
(pO| looks Uke 




Recall that in ([33| the constraints ^ > 0,b > 0, M > are supposed. The phase structure of the model in this 
particular case is defined by the properties of the global minimum point (GMP) of the TDP which certainly 

depend on the values of /i. The stationarity (gap) equations of this TDP, i.e. the equations dQ.^^'y^{M,b)/dM = 
and ai7Phy^(Af, b)/db = 0, read: 



Numerical investigations of the TDP (1551) and of the gap equations (|35p show that in the NJL2 model with 

chiral SUl{2) x SUii{2) symmetry the inhomogeneous CDW phase is realized only at /x > /Xc w 0.68 Mq. In contrast, 
at T = 0, in the (1+1 )-dimensional Ul{1) x Ur{1) chirally symmetric model the CDW phase appears at arbitrary 
nonzero values of /j, fVA [l8l | . Note that the critical value fic is equal to which corresponds to the lowest point of 
the homogeneous phase 2 (see Fig. 1). Below the critical chemical potential, i.e. at /Lt < fj,c, the usual homogeneous 
phase is arranged, where chiral symmetry is broken down to the diagonal SU (2) subgroup. The behavior of the chiral 
density wave amplitude M{fi) and its wave vector 60 (/^), which are the coordinates of the global minimum point of 
the TDP p3|) . is shown in Fig. 3 for /x/ — 0. It follows from this figure that at the critical point /ic a first order 
phase transition takes place, since here the order parameter M changes its value by a jump. Since in the CDW phase 
the relation M(^) < 60 (/x) < /x is valid, it is clear from the dispersion laws pT)) at A = that u-quarks are gapless 
excitations of this phase. It means that for each fi > fic there exist a momentum pi(/x) at which the quasiparticle 
energy pou is equal to zero, i.e. there is no energy cost to create u-quarks in the system. In contrast, for the energy 
of cZ-quarks we have throughout the CDW phase the relation pod > pomin = M{p) + boip) — ^ ~ M{fi), i.e. there is a 
gap in the energy spectrum of d-quarks which are called, due to this reason, gapped excitations of the CDW phase. 
There is one more peculiarity of the CDW phase. Indeed, as it is easily seen from ([33]) . at = the effective quark 
number chemical potential of u-quarks is equal to fi + b, whereas for d-quarks it is /i — b. Hence, there is a splitting 
of Fermi surfaces of up/down quarks by 2&o(m) in the CDW phase even at zero /i/. 

The fact that at fi = Mq the TDP ^ has a nontrivial minimum at the point (M w 0.25Mo, b w 0.99A/o) is well 
supported by figures Fig. 4, Fig. 5 and Fig. 6, where the plot of the function ^^^^^^{M, b) of M and b and its sections 
along axes b and M are presented (note that in order to draw the figures we continue the function p3l) symmetrically 
onto the negative semi-axis b). 

The influence of nonzero temperature on the formation of CDWs in the case fij — is considered, in particular, in 
the next section V (see Fig. 8). 



It is clear that to find the complete phase portrait of the model in terms of the external chemical potential 
parameters ly = ni/l and /i (at T = 0), one should investigate the global minimum point (GMP) of the physical TDP 
r2P'^y*'(M, 6, A) ((30)) vs the dynamical variables M,b,A. However, in the case under consideration the problem is 
simplified due to the effective reduction of external parameters. Indeed, the structure of ilP^y^{M, b, A) is such that 
it can be considered as a function of three dynamical variables Af, A,i> = b -\- v and only one external parameter /i, 
i.e. r2P''y''(Af, 6, A) = F{M, A, D; /i). So, the searching of the GMP of this function consists effectively of two stages. 
First, one can find the extremum of this function over M and A (taking into account the results of section III.B) 
and then, as it was done in the previous section IV. A, one minimizes the obtained expression over the variable v. 
Properties of the found GMP supply the following phase structure. 



As in the case with f) = 0, in the inhomogeneous case we did not found local minima of the TDP I I30I I of the form (M 7^ 0, A 7^ 0). 




(34) 



2b = 0{fi + b- M)y/{n + b)^ - AP +sign{b - n)e{\b M)^{b- ^if ^ NP. 



(35) 



B. General case: /ii ^ 0, /j. ^ 
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FIG. 7. The (i^, /i) phase portrait of the model at T = 



when a spatial CDW inhomogeneity is taken into account. 
In the CDWi (CDW2) phase 6 > (6 < 0). The curve L, 
on which & = 0, corresponds to the homogeneous chiral 
symmetry broken phase. The same is true for the interval 
< jU < /^a of the /i-axis, where /Xa = ~ O.68M0. 

If /X > /ict = /ic ~ O.68M0, then for arbitrary values of v we have found phases with spatially inhomogeneous 
condensates, which are realized at least in the form of chiral density waves or chiral spirals. The gap A is equal to 
zero for these phases. The amplitude M of these CDWs depends only on /i and is equal to the quantity M{n) (see 
Fig. 3). However, the chiral density wave vector b depends both on fi and v, namely 

b = bn{p) - V, (36) 

where the quantity 60 (/i) is also presented in Fig. 3. In the (:^, /i)-plane (see Fig. 7) we divide this region into two 
CDW phases. In the CDWi region we have the wave vector & > 0, i.e. here we have a clockwise twisted chiral spiral. 
In contrast, in the CDW2 region one obtains for chiral density waves the counterclockwise twisted chiral spirals, since 
here 6 < 0. For all points of the line L of this figure, which is defined by the relation L={(z/, ^) : v = 6o(/x)}, the 
wave vector b is equal to zero. So, the points of the curve L correspond to the homogeneous phase, where only chiral 
symmetry is spontaneously broken down and the dynamical quark mass is equal to M{fi) ^ (hence, on the line L the 
spatial translational invariancc of the system remains intact). Note, that the phase L is nothing else than the residue 
of the homogeneous phase 2 of Fig. 1 if the spatial inhomogeneity of chiral condensates is taken into account. To 
underline this fact, we use in Fig. 7 the notation fi^, which corresponds to the minimum point a of the homogeneous 
phase 2 of Fig. 1, for the critical curve between the CDW and charged pion condensation (PC) phases. However, fia 
coincides with the critical value of the case /x/ = (see section IV. A). 

As in the particular case with fij — (see the previous section), u-quarks are gapless excitations and d-quarks are 
gapped ones of the CDWi, 2 phases at fij ^ 0. The same is true for the homogeneous phase L. 

Below the line ji — oi Fig. 7 the homogeneous PC phase is arranged, since for all points of this region the GMP 
of the TDP (15(11) has the form M = 0, A = Mo,& = 0. In this phase the isospin Ui^{l) symmetry of the model is 
broken spontaneously. The exception is the interval < < /Xq, of the /i-axis, where chiral symmetry is broken down 
and quarks have the mass A/q. 

Note, both for the case of spatially homogeneous and inhomogeneous chiral condensate the isospin density n/ in 
the PC phase is equal to v/n. Starting from the VP^^'^^{M ,b, IS) it is possible to find the expression of this 

TDP in the CDWi, 2 phases (it is simply the expression (33) shifted by {—v^ /tt), in which Af, b should be replaced by 
M{^), boifJ.), correspondingly) and then to calculate their isospin density nj = —d^l^^^'^^/d^i. It turns out that in the 
CDW phases the isospin density is the same as in the PC phase, i.e. nj = v j-K. Hence, as it is easily seen from p6p . 
at fixed values of /x the wave vector of chiral spirals is tightly (linearly) connected with isospin density. In contrast, in 
the J7l(1) X Ur{X) symmetric NJL2 model without isospin chemical potential \ii the wave vector b shifts effectively 
the quark number chemical potential /x [13, [3 ■ Due to this reason, the quark number density Uq is equal to /x/tt in 
the CDW phase of this model. Moreover, the wave vector b in this phase is proportional to riq. 

V. CDW PHASES AT NONZERO TEMPERATURES 

In the case of spatially homogeneous condensates the infiuence of nonzero temperature on the phase structure of 
the S'?7l(2) X SUfj{2) symmetric NJL2 model (I) with two chemical potentials /i and v = /x//2 was considered in |24| . 
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Now let us study the influence of temperature T on the phase structure of this model in the case of an inhomogeneous 
chiral condensate of the form In this case, to get the corresponding (unrenormalized) thermodynamic potential 
Q.T {M, b, A) one can simply start from the expression for the TDP at zero temperature (|15|) and perform the following 
standard replacements: 



PO POn 



i7rT(2n + l), n = 0, ±1, ±2, 



(37) 



n— — oo 



i.e. the po-integration should be replaced by the summation over an infinite set of Matsubara frequencies uj„. Summing 
over Matsubara frequencies in the obtained expression (the corresponding technique is presented, e.g., in [ssj), one 
can find for the TDP: 



Ot(M,6, A): 



AP + A^ 
4G 

Tin [1 + 



r ^|i5+ + ij^+Tln[l + e-^(^i-^)] +Tln[l + e-^(^i+'')] 



-fi{,E^-^J.)^ 



Tin [1 



(38) 



where P = l/T and are given in ([TB)) . Clearly, only the first two terms (which are the same as in the zero 
temperature case), in the braces of this expression are responsible for an ultraviolet divergency of the whole TDP 
([55)1 . So, regularizing the TDP ([55]) in the way as it was done in (l5lT) for zero temperature TDP and then replacing 
G — !■ G(A) (see formula (j20p ). we can obtain in the limit A — > oo a finite expression denoted as 17f^"(M,6,A). It IS an 
evident generahzation of the TDP fiP'^y"(M, 6, A) ^ to the case of nonzero temperature. Numerical investigations 
show that all possible local minima of the obtained TDP rjP''y'(M, h, A) are located in the planes M = or A = 0. 
So it is sufficient to deal with corresponding restrictions of the TDP on these planes, i.e. with the following functions. 



nf^'^iM - 0,6, A 
r!fy^(M,6,A 



V^,(A) - 
0)=Fo(M) 



2T 

TT 



dpi In • 



-/3(£-m) 



by 



T 



dpi In I 



1 + e 



-/3(£+m) 



T 

TT 



dpi In I 



1 + e 



1 + e 



1 + e-/9(-E+'v+6+M) 



(39) 



(40) 



where the effective potential Vq{x) is given in (|2ip . E ~ ^/j^+M^, and £ = ^/pf + A^ . Comparing the global 
minima of the functions (j39]) and (|40|) , it is possible to establish the global minimum point of the renormalized TDP 
fl^^^^{M, b, A). Then, the dependence of the global minimum point vs T, /i, v defines the phase structure of the model. 

Using this prescription in our numerical investigations of the TDPs ([55)1 - (|in)) . we have found the two (/i,T)-phase 
portraits of the initial NJL2 model (1) depicted in Figs. 8, 9 for qualitatively different fixed values of the isospin 
chemical potentials, < /i/ < 2j/q and 2i^q < A*/, respectively ~ O.6M0 is the i^-coordinate of the point a of 
Fig. 1). Note, there is a first order phase transition on the boundaries between CDW1.2 and homogeneous PC or 
chiral symmetry breaking phases of these figures. However, other boundaries of the phases of Figs. 8, 9 correspond 
to critical curves of the second order. It is interesting to remark that for < v < i^a {^a < v) all critical curves of 
Fig. 8 (Fig. 9) do not depend on v. 

Finally, let us take fj,i = and compare the thermodynamical properties of our (l+l)-dimensional NJL model 
(1) (see the phase portrait of Fig. 8 at = 0) with the corresponding massless (3+l)-dimensional NJL model with 
chiral S'[/l(2) x SUr{2) symmetry [29|. It turns out that in the four-dimensional spacetime, in contrast to the (1+1)- 
dimensional case, a second order phase transition from a homogeneous chirally broken phase to an inhomogeneous 
one takes place. Moreover, depending on the value of the dynamical quark mass in vacuum, the inhomogeneous 
phase in the (3+l)-dimensional NJL model may occupy both a finite (compact) and infinite (noncompact) region 
of the (//, r)-phase diagram, whereas in our two-dimensional NJL model (1) an inhomogeneous phase appears as a 
noncompact region (see Fig. 8). 



VI. CONCLUSIONS 



It is well known that at nonzero baryon densities there might exist phases with a spatially inhomogeneous chiral 
condensate. This fact was established in the framework of both two-dimensional GN- or NJL-type models [T7l.[l8l[26j 
and four-dimensional NJL-type models [27l - [3H . where phases with a crystalline chiral condensate or with a CDW 
spatial inhomogeneity were proved to exist at nonzero values of the baryon chemical potential. Since the isotopic 
asymmetry of dense quark matter is an inevitable reality, a more adequate investigation of dense quark matter 
demands to include into consideration both baryon, ^, and isospin, /i/, chemical potentials. In this paper and in 
contrast to previous papers [13, [H, , we study the possibility of spatially inhomogeneous chiral condensates in the 
SUl{2) X SUr{2) symmetric NJL2 model (1) including the two above-mentioned chemical potentials in the leading 
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FIG. 8. The {n, T) phase portrait of the model at fixed 
V, where Q < v < Va O.GMq. There, at = H denotes 
the homogeneous chiral symmetry breaking phase with 
M = Mo,b = 0, A = 0. At < ly < i/a U denotes 
the homogeneous charged pion condensation phase (PC) 
with M = 0, 6 = 0, A = Mo. fic ^ fJ-a ^ O.68M0. In the 
symmetric phase M = 0, b = 0, A = 0. CDWi denotes 
an inhomogeneous chiral density wave phase with fe > 0. 
All critical curves do not depend on v. 



T/Mo 



SYMMETKIC PHASE 




FIG. 9. The (/i, T) phase portrait at fixed v, where 
Va < I'- Here bQ^{i') is the function inverse to &o(m) 
defined in Fig. 3. PC denotes the homogeneous charged 
pion condensation phase with AI — 0,b = 0,A = Mq. 
CDW2 denotes the inhomogeneous chiral density wave 
phase with 6 < 0. All critical curves do not depend on 1/. 
Other notations are the same as in the previous figures. 



order of the large- TVc expansion. For simplicity, the spatial inhomogeneity in our consideration is realized in the form 
of so-called chiral density waves or chiral spirals. 

First, we have proven that at — and T = the inhomogeneous CDW phase is realized in this SUl{2) x SUr{2) 
symmetric NJL2 model only at sufficiently large values of ^, i.e. dX jjL > ~ O.68M0 (here Mq is the dynamical 
quark mass in the vacuum, i.e. at zero values of chemical potentials). In contrast, it is well-known that in the NJL2 
model with continuous ?7l(1) x Ur{1) chiral symmetry the CDW phase appears at arbitrary nonzero values of /i > 
[l3, Moreover, it turns out that at /z/ = the Fermi surfaces of up/down quarks in the CDW phase are split by 
26o(m)) where 60 (/^) is the wave vector in this phase. 

Second, if /i/ ^ and T = then in the (/i/,/i) phase diagram (see Fig. 7) the spatially inhomogeneous phases 
are allowed at fi > fia — fJ-c and arbitrary values of fij. This region is divided by the curve L into two domains. 
In one of them each CDW is a clockwise twisted chiral spiral, in the other - it is a counterclockwise twisted spiral. 
The amplitude of chiral density waves does not depend on fij. The dependence of its wave vector b on fi and fij is 
defined by the formula (1361) . Since the isospin density nj in these phases is equal to v/ir, we see that the wave vector 
b is linearly connected with nj. In contrast, in the Uijl) x [/ii;(l)-symmetric NJL2 model the wave vector of the 
CDW phase is proportional to a quark number density [13, [3 • The points of the curve L correspond to the spatially 
homogeneous phase (since here 6 = 0) with spontaneous chiral symmetry breaking. Indeed, the phase L is the residue 
of the homogeneous massive chirally nonsymmetric phase 2 of Fig. 1 which shrinks to L after taking into account 
inhomogeneity phenomena. Below the line fi = fia the homogeneous charged pion condensation phase is realized. 

It turns out that at arbitrary /^/-values in all above mentioned inhomogeneous CDW phases as well as in the L 
phase li-quarks are gapless excitations, but d-quarks are gapped ones. 

Third, we have studied the influence of temperature on the formation of the CDW phases. In particular, it was 
shown that at /i/ = the (/i,T)-phase diagrams of the SUl{'2) x SUii{2)- and C/l(1) x Ur{1) symmetric NJL2 models 
are quite different. Indeed, as it was proved in Ts], in the second model the CDW phase occupies in this diagram an 
infinite strip which includes points with arbitrary small /i-values, whereas in the first model (see Fig. 8) the upper 
boundary of this phase is a monotonically decreasing function of /i. In addition, for rather small values of ^ the CDW 
phase is forbidden in the framework of the SUl{2) x SUr{2) symmetric NJL2 model. 

We finally note that in this paper we have suggested a homogeneous pion condensate. It would be interesting to 
study in future the possibility of the spatially inhomogeneous pion condensation phase. 
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